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DISCRETE TOMOGRAPHY SOLUTIONS 



LAJOS HAJDU AND ROB TIJDEMAN 

Abstract. In earlier papers we have developed an algebraic the- 
ory of discrete tomography. In those papers the structure of the 
functions / : A — > {0, 1} and / : A — > Z having given line sums 
in certain directions have been analyzed. Here A was a block in 
Z" with sides parallel to the axes. In the present paper we assume 
that there is noise in the measurements and (only) that A is an 
arbitrary or convex finite set in Z n . We derive generalizations of 
earlier results. Furthermore we apply a method of Beck and Fiala 
to obtain results of he following type: if the line sums in k direc- 
tions of a function h : A — > [0,1] are known, then there exists a 
function / : A —> {0, 1} such that its line sums differ by at most k 
from the corresponding line sums of h. 



1. Introduction 

Let n be a positive integer and let A be a finite subset of Z n . If 
/ : A — > R, then the line sum of / along the line / = c + td (with 
c, d G Z n , d 7^ fixed and £ G IR variable) is defined as J2 a eAni fiod- 
We call d a direction. Let S = {di, . . . , d^} be a set of directions. By 
the line sums of / along S we mean all the line sums of / along a line 
in a direction from S passing through at least one point of A. Theorem 
1 of [H] states that if A is a block with sides parallel to the axes, then 
any function / : A — > K. with zero line sums along S can be uniquely 
written as a linear combination of so-called switching components of S 
contained in A. In Section 3 we prove that for this result it suffices that 
A is convex, but that the convexity requirement cannot be dropped. 

By a discrete tomography problem we mean asking for a function 
/ : A — >■ Z which satisfies prescribed line sums along S, where Z may 
be {0,1}, M, Z or some finite real set. The authors and others have 
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developed an algebraic theory of the structure of the solutions of a 
discrete tomography problem, see [8], [9], [ID], [EE], 0) [13]) [13]) 0) 
[3]. It appears that the real solutions of a discrete tomography problem 
form a linear manifold if there is at least one real solution, and that 
the integer solutions form a grid in this linear manifold, provided that 
at least one integer solution exists. 

If the line sums are measured with some noise, then it is not certain 
that some function satisfies the measured line sums along S. A natural 
question is then what the best approximative solution is. We shall show 
that there is some linear manifold which can be considered as the set 
of 'best real approximations' in the sense of least squares. An obvious 
choice is then to choose the shortest best approximation, that is the 
orthogonal projection of the origin to that linear manifold. In Section 4 
we present an algorithm to construct this shortest best approximation 
and illustrate it by an example. In Section 5 we present an explicit 
system of linear equations which determines the shortest best solution 
in case A is convex. As an application we generalize a result from [5] 
by giving an explicit expression for the shortest best solution in case A 
is a rectangle with sides parallel to the axes and only row and column 
sums are given. 

In the 80's Beck and Fiala [I] proved a 'balancing' theorem. In 
Section 6 we show that this implies that if the line sums in k directions 
of a function h : A — > [0, 1] are known, then there exists a function 
/ : A — > {0, 1} such that its line sums differ by at most k from the 
corresponding line sums of h. 

We extend this result in Section 7 to the case that we are not search- 
ing for a binary image, but for an image / with a finite number of given 
real values. To do so we generalize the result of Beck and Fiala. 



2. Notation 

We use the following notation throughout the paper. Let n be a 
positive integer. For brevity, for x\, . . . , x n G R and u\, . . . , u n G Z, we 
write x = (xi, . . . , x n ), x = (xi, . . . , x n ) T and x- = YYj=i x / ■ 

Let d G Z n with gcd(c?i, . . . , d n ) = 1 be such that d ^ 0, and for the 
smallest j with dj ^ we have dj > 0. We call d a direction. By lines 
with direction d we mean lines of the form c + td (with c G II 1 fixed, 
t G R variable). 

Let A be a finite subset of Z n . Write A = {ai,...,a s } where 
ai,...,a s are arranged in lexicographic increasing order. We call A 
convex if every a G Z n which belongs to the closed convex hull of A 
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belongs to A itself. By the minimal corner of a set B C Z n we mean 
the lexicographically smallest element (f>(B) of B. 

If / : A — > M, then the line sum of / along the line / = c+td is defined 

as EaeAnz /(e)- For an y / : A -> R ' write / := (/(fii)> • • • > f(a s )) T - 
We often identify / and /. The length of / (or /) is defined as |/| = 

1/1 = y/Ea^im) 2 - 

Let k be a positive integer and 5 = {d\, . . . , dfc} be a fixed set of 
directions. By the line sums along S we mean all the line sums along 
lines in a direction from S which pass through at least one point of A. 
For d = (di, . . . ,d n ) G S put 

^) = fe--l)II^. 

dj<0 

= Yli=i fdA%) an d, for u G Z n , set FJ^x) = x-F(x). Obviously, 
the polynomial has integer coefficients. We call the functions F^ 
the switching polynomials of S. Define the functions m„ : Z n — > Z by 

m>u{v) = coeff(x-) in for v G Z n . 

We define D„ as the set of v G Z™ for which m„(t)) 7^ and call it a 
switching component. Let 0(V) denote the minimal corner of D„. It 
follows from the above definitions that 

(1) mu{(j>{u)) = ±1. 

3. The structure of functions with zero line sums 

We prove that Theorem 1 of [TTJ remains true under the weaker 
condition that A is convex. 

Theorem 3.1. Let A be a finite convex subset of Z n , and S a given 
set of directions. Then any function f : A — > R with zero line sums 
along S can be uniquely written in the form 

DuCA 

with coefficients Cu G M. Moreover, every such function f has zero line 
sums along S . 

If there is no u for which L>„ C A, then the only function / with zero 
line sums along S is the trivial function / = 0. Otherwise, the functions 
with zero line sums along S form a proper linear subspace of the linear 
space of all functions / : A — Y R. 
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Proof. The statement has been proved in case A is a hyperblock with 
sides parallel to the axes in Theorem 1 of [UJ . Let A* be a hyperblock 
with sides parallel to the axes such that A C A*. Set f(x) = for 
x G A* \ A. Then we know that 

(2) / = c„m„ 

D„CA* 

with coefficients c u G R. It remains to prove that c u = if _D U is not 
contained in A. 

If -D„ is not contained in A, then there exists ip(u) G .D^ such that 
Tp(u) ^4- Since A is convex, there is a linear manifold L which extends 
a hyperface of the convex hull of A such that ip{u) and A are on different 
sides of L. Let Hl be the open halfspace generated by L which contains 
i/j(u). Note that Hl does not contain any element of A. Consider the set 
Ul of all u such that Du Q A* and Du contains an element if)(u) G Hl- 
Without loss of generality we assume that ip(u) has maximal Euclidean 
distance d(ijj(u), L) to L among the elements of Du D Hl and, if there 
are more such elements with maximal distance to L, then ip(u) is the 
lexicographically smallest among them. Since the sets Du for variable 
u are translates of each other, the vectors ip(u) — u are the same for 
all u(zUl. Now we arrange the elements of Ul according to the non- 
increasing distances d(ip(u),L) of ip{u) to L. Thereafter we order the 
elements of Ul for which the distances d(ip(u), L) are equal according 
to non-decreasing lexicographic order of u. Consider the first element 
u&Ul according to this ordering. By the above construction there is 
no other set D„ for u G Ul which contains i/j(u). Since ip(u) A we 
infer f(ip(u)) = 0, hence c M = 0. We proceed with the next element 
u G Ul in the ordering and conclude by a similar reasoning that = 
also for this u. Continuing until we have had all elements of Ul-, we 
conclude that c„ = for all u^Ul- Since D„ was an arbitrary set not 
contained in A, the first statement follows. The uniqueness and the 
second statement of the theorem follow immediately from Theorem 1 
of P3]. □ 

The following result is a consequence of Theorem 13.11 

Corollary 3.1. In the notation of Theorem \3.1\ for any h : A — > R and 

for any prescribed values from R at the minimal corners of the switching 
components contained in A there exists a unique f : A — > R having the 
same line sums along S as h has and having the prescribed values at 
the minimal corners. Moreover, if h : A — > Z, then f : A — >■ Z. 
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Proof. According to Theorem 13.11 there are unique coefficients c„ such 
that 

/ = h + 22 c « m « 

has the same line sums along S as h. By ([1]) we obtain, following the 
ordering argument from the previous proof, that each coefficient c„ is 
completely determined by the value of m„ at (f>(u) and, moreover, that 
Cu e Z if h : A -> Z. ~ □ 

Remark 3.1. The following example shows that in Theorem 3.1 we 
cannot drop the convexity requirement. 

Let A = {(0, 0), (0, 1), (1, 0), (1, 2), (2, 1), (2, 2)} and S = {(1, 0), (0, 1)}. 
Then for every u we have D % — u = {(0, 0), (0, 1), (1, 0), (1, 1)}. There- 
fore A does not contain any switchting component. However, there is 
a nontrivial function / : A — > Z with all line sums along S equal to 0: 
/(0, 0) = 1, /(0, 1) = -1, /(l, 0) = -1, /(l, 2) = 1, f(2, 1) = 1, 
/(2,2) = -l. 

4. The best approximating function for general domains 

The next theorem can be used to construct the function f : A — > K. 
such that /o fits optimally the measured line sums along 5 in the sense 
of least squares and, moreover, has minimal Euclidean length among 
such functions. 

Let A C Z n be a finite, nonempty set, and write a x , . . . ,a s for its 
elements. For the rest of this section, fix the indexing of the elements. 

Let B be a t by s matrix of real numbers. The range of the matrix 
B is denoted by 

R{B) :={B-x : x G 1R S }. 
Hence -R(-B) is a subspace of R*, generated by the column vectors 
6 a of 5. We have < dim( J R( J B)) < £. Write B x for a ma- 
trix formed by a maximal linearly independent set of column vectors 
of B. Then B\ = B ■ C\ where C\ is a matrix of type s x (rank(S)) 
which has rank(5) entries 1 in distinct columns and all other entries 
equal to 0. Observe that Bf ■ B\ is invertible. 

Lemma 4.1. Let A,B and B\ be as above. Let h e M* be arbitrary. 
Put 

(3) b* = B x ■ {B\ ■ B,)- 1 ■ Bl ■ b. 

Then b* is the vector from R{B) which is closest to b. 

Proof. Obviously, the vector b* in R(B) closest to b is uniquely deter- 
mined by the following properties: 
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• 6* e R(B), 

• b — b* is orthogonal to R(B). 

Since B\ = B ■ C\ the first property follows immediately from fl3]) . The 
second property is equivalent to that b — b* is orthogonal to all the 
column vectors of B, or equivalently, to all column vectors of B\. In 
other words, it is equivalent to 

B{ ■ (6-6*) = 0. 

It follows from fl3]) that 

Bl ■b* = Bl-B l - (Bj ■ BJ- 1 -Bl-b = Bl-b. 
Hence both properties are satisfied. □ 
We use the above notation and define 

h = B ■ f. 

Let B2 be a matrix formed by a maximal linearly independent set of 
row vectors of B. Then B2 = C2 • B where C2 is a matrix of type 
(rank(I?)) x t which has rank(5) entries 1 in distinct rows and all 
other entries equal to 0. Observe that B 2 ■ B% is invertible. 

Theorem 4.1. Let A, B, B 2 , C 2 , 6, 6*, / (= f) be as above. Put 

(4) f = Bl-{B 2 -BlY l -C2-b* 

Then the corresponding f : A — >■ K has the following properties: 

(i) for any f : A — > R we have \lf — 6| > \lf — b\, 

(ii) if f : A — > R, / 7^ / suc/i f/iai |// — b \ = \lf — b\, then 
\f\>\fo\. 

Proof. Observe that (i) and (ii) are equivalent with the following two 
properties: 

• B-fo = b*, 

• fo is orthogonal to ker(B), the nullspace of B. 
The first property is clearly equivalent to 

(5) B 2 -f = C 2 - b*. 

It follows immediately from that this property is satisfied. Since 

ker(B) = {xeR s : B-x = 6}, 

the orthogonal complement of ker(B) is the subspace of R s generated 
by the row vectors of B. Hence, by the definition of B 2 , the second 
property above is equivalent to 

(6) fo = B2 ■ y for some y G R r , 
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where r is the rank of B 2 - This is obviously true because of (BJ. Thus 
both properties are satisfied. □ 

Remark 4.1. An alternative version of Theorem 14.11 can be obtained 
by using the Moore-Penrose pseudo inverse, cf. the proof of Theorem 
1 in [3]. 

Remark 4.2. We apply Lemma 147X1 and Theorem 14. II in the context of 
Discrete Tomography as follows. Let A be a finite subset of Z n and S a 
set of directions. Let l lf . . . , l t be the measured line sums along S. Note 
that because of noise they need not be consistent. Then B is the s by 
t matrix whose entry Bij equals 1 if the line corresponding to lj passes 
through etj and otherwise. The vector b* constructed in Lemma 14.11 
represents the corresponding line sums along S which are consistent and 
provide the optimal choice in the sense that Y^j=i(^j ~ tf) 2 * s m i mm al 
among the consistent line sums b* along 5*. Furthermore, the vector 

fo constructed in Theorem 14.11 is the shortest best approximation in 
the sense that it is the shortest vector realizing the line sums given by 
b*. The corresponding function f : A — >■ R may be considered as the 
optimal choice for the measured line sums l±, . . . ,l t . 

We illustrate the method by an example. 

Example. We use the notation of Lemma [4.11 and Theorem 14.11 Con- 
sider the following subset of M 2 : 

A := {(1, 0), (3, 0), (0, 1), (4, 1), (0, 2), (4, 2), (1, 3), (2, 3), (3, 3)}. 

As the set of directions, take 

S:= {(1,0), (0,1), (1,-1), (1,1)}. 

The ordering of the points in A and directions in S are arbitrary, but 
fixed. As a (measured) line sum vector, take 

' V 10' 5' ' ' '2' '5' ' ' '10' 10' 2' '5'5'2' 10' 10; ' 
The entries of b belong to the lines 

y = t (t = 0,1,2,3), x = t (t = 0,1,2,3,4) 



y = x + t (£ = -3,-2,-1,0,1,2), y = -x + t (t = 1, 2, 3, 4, 5, 6) 
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which we keep in this order. Then the matrix B of line sums is given 
by 
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As one can easily check, rank(i?) = 9. So we can take the matrix C\ 

as the 9x9 unit matrix. Thus Bi = B. Then, by ([3]), the vector b* 
is given by 



I 891 2457 1019 4361 167 103 111 103 963 4239 859 1211 1433 287 2511 4239 1179 571 153 367 3151 \ 
I 800 ' 1600 ' 800 ' 3200 ' 128 ' 128 ' 128 ' 128 ' 640 ' 3200 ' 1600 ' 1600 ' 3200 ' 200 ' 3200 ' 3200 ' 1600 ' 1600 ' 3200 ' 200 ' 3200 ) 1 

As one can readily check, the indices of a maximal set of independent 
rows of B is given by 

{1,2,3,4,5,6,7,10,11}. 

That is, we may take 
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whence 



Bo 



/1 1 0000000\ 

1 1 

1 1 

1 1 1 

1 1 

1 1 








1 




1 1 
\000001000/ 

Finally, by (jl]) we obtain 

/1211 571 1817 3097 1179 859 153 111 1433 \ 

V1600' 1600' 3200' 3200' 1600' 1600' 3200' 128' 3200 J 



fo 



5. The best approximating function for convex domains 

The following theorem provides explicitly a system of linear equa- 
tions which determines the best approximating function constructed in 
the previous section. We illustrate in the corollary the advantage of 
this explicit expression. The real number l(Y T ) in the following the- 
orem can be considered as the measured line sum of / along the line 
corresponding to Y T . 

Theorem 5.1. Let A C M ra be convex. Let S be a finite set of direc- 
tions and Y\, . . . ,Y t the subsets of A which determine the lines along 
S . Suppose for r = 1, . . . , t a real number l(Y T ) is given. Let Ua C A 
be the set of minimal corners of the switching components contained in 
A. Define f :A^-M,by the system of linear equations 



(7) 



o(v) m u(v) = for all u with G Ua, 



(8) £ f°(^ = l ( Y ^ f° r aU & mth ^) G A \ Ua - 

T:U£Y t V_& T T:ueY T 

Then f is a function such that 



(9) £ £/ofe)-i(n) 

T = l \veY T 

is minimal and among such functions fo is the one for which the value 
°f I fo I is minimal. 

Proof. By Theorem 14. li the function fo : A — > R satisfying for which 

\fo\ 2 = J2veA(fo(^-)) 2 ^ s minimal is uniquely determined. We proceed 
with this function / and consider it as a function for which each value 
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fo(u) for u G A is a variable. It follows by different at ion of (jUJ) to each 
f (u) that 

E E /ooo = E to 

t:u6"Kt ueK T t:u£Y t 

for all uE A. Hence / satisfies (E}. 

We know that fo is orthogonal to the linear subspace L of functions 
having zero line sums along S. According to Theorem 13. II the functions 
rriu have zero line sums along S. Therefore they are in L for all aeZ" 
Since the inner product of fo and any vector from L is 0, fo satisfies 
© too. 

The numbers of linear equations in ([7]) and flB]) together equal the 
cardinality of A. Thus it suffices to show that they are linearly inde- 
pendent over R in order to prove that fo is completely determined by 
them. Because of the orthogonality of f and L, it is enough to prove 
that the equations in ([7]) are linearly independent as well as those in 
®. 

Since by Theorem 13. II the functions m„ are linearly independent, the 
equations ([7]) are linearly independent as well. 

Furthermore, in Theorem 13.11 it is shown that fo is uniquely deter- 
mined by its values at Ua- This shows that the equations in (JE]) are 
linearly independent. We conclude that the linear equations in ([7]) and 
([HD are linearly independent indeed. □ 

In the particular case that A C Z 2 is a rectangular block, and we 
only have row and column sums, we give an explicit form of fo- The 
result shows that the formula from [6J is also valid if there is noise in 
the measurements. We simplify our notation. 

Corollary 5.1. Let A = G Z 2 : < i < q, < j < p}, S = 

{(1, 0), (0, 1)}. Let Ci (i — 0, . . . , q — 1) and Tj (j — 0, . . . , p — 1) denote 
the measured column sums and row sums, respectively. Further, write 

p-l q-l 

s r = r ji s c — S c i an d T = pSr q +p Sc ■ Then for any {i,j) G A we 



j=0 i=Q 

have 



f (i,j) = - + r - 1 --. 

p q qp 



Observe that if s r = s c , then T = s r = s c . 
Proof. Since 

{ J -+- — - -+^ ±i — )-(—+- — + M ±i +^ ±i — 

q p qp q p qp q p qp q p qp 
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for all i and j, the equations (EJ) are satisfied. Furthermore 

( a+£i_l )+ ... +( i> + Si_l )+( ii + ^_l )+ ... +( a + 5s_l ) 

q p qp q p qp q p qp q p qp 

= h Q h Tj H = Cj + Tj, 

q p VI 

which shows that the equations (jSD are also satisfied. □ 

6. Approximate solutions in the binary case 

Let A be a finite subset of Z n . We assume that a function h : A — > M 
is given and provide information on the 'nearest' function / : A — > Z 
having approximately the same line sums along S as h. 

If n = 2 and only row and column sums are given, we have the 
following result. 

Theorem 6.1. Ifh : A — > R is given, there exists a function f : A — >■ Z 
s«c/i i/iai ever?/ two corresponding elements of f and h as well as every 
two corresponding row sums and column sums as well as the sums of 
all function values of f and h differ by less than 1 . 

We apply the following result of Baranyai. 

Lemma 6.1 ([2], Lemma 3). Let [hij] be an I by m matrix of real 
elements. Then there exists an I by m integer matrix [f\j] such that 

\hij - fij\ < 1 for all 



i i 

I hij ~ ^2 fa I < 1 for a11 

3 j 
t 3 i j 

Proof of Theorem \6.1[ Choose an I by m block A* which covers A. For 
G A* \ A put h(i,j) = 0. This does not change the line sums. 
Applying Lemma IfTTT we get f(i,j) = h(i,j) = for G A* \ A 
and the theorem follows. □ 

The following example shows that the bound 1 is best possible. Let 

< s < 1, I > 1/s, m = 1, h(i, 1) = e for i = 1, . . . , I. Then f(i, 1) = 1 
for some i in order to avoid that the row sums of h and / differ more 
than 1. But then the z-th column sums of h and / have a difference 

1 -£. 
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The crucial feature of the following general result is that the upper 
bound is independent of the size of A. 

Theorem 6.2. Let A be a finite set in Z n . Let h : A — > R and let k 

directions S be given. Then there exists a function f : A — )■ Z such 
that each difference between corresponding elements of h and f is less 
than 1 and each difference between corresponding line sums of h and f 
along S is at most k — 1. 

We introduce the following notation in order to apply a result of 
Beck and Fiala. Let X = {x±, X2, ■ ■ ■} be a finite set and J 7 a family 
of subsets of X. Associate to every Xi a real number c^. Let k be the 
degree of J 7 , that is the maximal number of elements of T to which 
some element of X belongs. Let r(k) be the least value for which one 
can find integers a iy i = 1, 2, . . . so that |a» — a»| < 1 and 

for all E G J 7 . The following result is due to Beck and Fiala (see [1]). 
We shall prove a generalization of it in the next section. 

Lemma 6.2. In the above notation, we have 

r(k) < k - 1 for k > 2. 

Beck and Fiala conjecture that r(k) < k/2 is true even for small values 
of k. Bednarchak and Helm [5] and Helm [T2] improved the Beck-Fiala 
bound to r(k) < k — 3/2 for k > 3 and r(k) < k — 2 for k sufficiently 
large, respectively. 

Proof of Theorem W.B. Let Yi, . . . ,Y t denote the subsets of A which de- 
termine the line sums along S. Let J 7 = {Y\, Y2, . . . , Y t }. By Lemma 
16.21 there exist integers f(a) for all a G A with f(a) G {|_M a )_K TM a )l} 
such that Ylia&j \f( a ) ~ h(a)\ < k — 1 for j = 1, . . . , t. □ 

Remark 6.1. Obviously, many variations of Theorem 16 .21 are possible. 
E.g. adding the requirement that the sum of all values f(a) differs 
little from the sum of all values h(a) leads to an upper bound k. The 
requirement that the difference between the sums of the values of / 
and h along any linear manifold parallel to the axes should be small 
leads to an upper bound 2 k — 2. 

Remark 6.2. By a probabilistic method a better dependence on k can 
be obtained at the cost of some dependence on A. An recent improve- 
ment by Banaszczyk [1] of a result of Beck implies that in Theorem 16.21 
the upper bound k — 1 can be replaced by C A/Hog(min(m, n)), where 
C is some constant. 
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7. Approximate solutions for grey values 

Theorem 7.1. Let Z = {zi, . . . , z m } be a set of m real numbers with 
Zi < • • • < z m . Put z = maxi(z i+ i — zA. Let h : A — >■ [z±, z m ] and let 
k directions S be given. Then there exists a function f : A — >■ Z such 
that the difference between the values of f and h at any element of A 
is at most z and each difference between corresponding line sums of f 
and h along S is at most (k — l)z. 

For the proof we derive the following extension of the lemma of Beck 
and Fiala. 

Lemma 7.1. Let Z = {z\, . . . , z m } be a set of m real numbers with 
z± < • • • < z m . Put z = m&Xi(z i+ i — z^. Let X = {xi,x 2 , ■ ■ ■ ,x s } be 
a finite set and associate to every Xi a real number cti G [zi, z m ]. Then 
given any family T of subsets of X having maximum degree k > 2, 
there exist a\ G Z such that ai = Zj if a, 
from Z in between cti and for all i and j and 



Zj and there is no element 



< (k - l)z 



for all E G T . 



Proof. We shall define a sequence a , a 1 , . . . , a p of s-dimensional vec- 
tors ai = (a{, . . . , a J s ) and a sequence Yj of subsets of X with the 
following properties: 

i) = ctj for i — 1, . . . , s. 

ii) There is no element of Z in between ctj and a\ for % — 1, . . . , s; j — 
>,!,•■ -,P- 

iii) X \ Yj is a set of points x for which x G Z for all j. 

iv) Y D Yi D • • • D Y p and \Yj\ = p - j for < j < p. 

v) a\ = a\ for j = h, . . . , p whenever a\ G Z . 

vi) If \EC\Yj\ > k, then J2 X 



EE 



a:-, 



E, ieE ^' +1 for all E E 



vii) For j — 0, 1, . . . , p and all E G J 7 we have 



a. 



Xi£E 



XieE 



< (k - l)z. 



According to (iii) and (iv) the final vector a p has all coordinates in Z. 

We construct the sequence (a- 7 ) by induction. Suppose is defined 
satisfying the above conditions for j. Let 

Gj = {E G 7 : \ECiYj\ > k}. 
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We distinguish between three cases. At every step there is some i such 

that Xi G Yj, a{ +1 G Z and we set Yj +1 = Yj \ {xi}. 

Case (a) Gj = 0. 

Case (b) < \Gj\ < \Yj\. 

Case (c) \Gj\ > \Yj\. 

Case (a). If Gj is empty, then choose a{ +1 as the element from Z 
which is nearest to on for all i with x^ G Yj. It follows that 



oa - J2 a i +1 



<(k — l)z for all E G J 7 , 



and the above conditions are satisfied for j + 1. 
(It follows that a{ +l — ■ ■ ■ = of = £Jj for all z.) 

In Case (b) associate a real variable to every z = l,...,s and 
consider the system of equations 

^ A = for E G Gj, 

(3i = for zrj ^ Y}. 

A nontrivial solution {/3i}f =1 exists, because in case (b) there are more 
variables than equations. Let to be the smallest nonnegative value for 
which a{ + tfli G Z for some % with Xi G Yj. Put a{ +1 = a\ + to A for 
i = 1, . . . , s. It is easy to check that 

^2 <4 = J2 a i +1 for al1 EeGj. 

Xi£E Xi£E 

Hence the above conditions are satisfied for j + 1. 

Case (c). Since each Xi has degree at most k in Gj, we may conclude 
that \Gj\ = \Yj\, each Xj has degree exactly k in Gj and \E n V}| = 
for every G Gj. Let a{ +1 be the element from Z nearest to aj for 
every Xj G Y/. Then \a{ +1 — ati\ < z/2 for Xi G Y}. Since k/2 < k — 1, 
we obtain 

|^«r-^« 4 |<(fc-l) Z 

XiGE Xi£E 

for all G J 7 . Hence the above conditions are satisfied for j + 1. 
(It follows that a{ +1 — ■ ■ ■ — a? — a* for all i.) 

Write <jj = a? for i = 1, . . . , s. It is easy to check that in each case 
the relations (iii), (iv) and (vii) hold. This completes the proof. □ 

Proof of Theorem \ 7. 1\ Let Y"i, . . . , Yj denote the subsets of A which de- 
termine the line sums along S. By Lemma 17.11 there exists a function 
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/ : A — > Z such that 

\f(a)-h(a)\<(k-l)z 

aGAnYj 

for j = l,...,t. □ 

Remark 7.1. A small adjustment must be made if the entries are not 
all in [zi,z m ]. E.g. values of h smaller than z\ are first replaced by z\, 
values larger than z m by z m . 

Remark 7.2. If we want to have relatively short vectors f,g, then we 
may apply Theorem 17.11 to the function fo from Theorem 14. 1[ 
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